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Abstract 

We study the full counting statistics of interferometers for chiral Majorana 
fermions with two incoming and two outgoing Dirac fermion channels. In 
the absence of interactions, the FCS can be obtained from the 4x4 scatter¬ 
ing matrix S that relates the outgoing Dirac fermions to the incoming Dirac 
fermions. After presenting explicit expressions for the higher-order current 
correlations for a modihed Banbury Brown-Twiss interferometer, we note 
that the cumulant-generating function can be interpreted such that unit- 
charge transfer processes correspond to two independent half-charge transfer 
processes, or alternatively, to two independent electron-hole conversion pro¬ 
cesses. By a combination of analytical and numerical approaches, we verify 
that this factorization property holds for a general scattering matrix, 

i.e. for a general interferometer geometry. 
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1. Introduction 


In a seminal work Biittiker pointed out the power of investigating non¬ 
local current correlations in mesoscopic conductors to detect particle ex¬ 
change effects [T]. This work extended the scope of previous studies of shot 
noise in two-terminal conductors [SlElilE] by showing that fundamental 
quantum mechanical effects can play a decisive role in electronic transport 
properties. 

Later on, Biittiker demonstrated on completely general grounds that 
cross-correlations of bosonic and fermionic free particles are fundamentally 
different |6]. Whereas for bosons the sign is determined by the competition 
between positively contributing bunching effects and negatively contribut¬ 
ing partitioning effects, free fermions always exhibit overall negative cross¬ 
correlations due to their antibunching property. Thus, the observation of 
positive cross-correlations for fermions requires some nontrivial interaction 
between them. A hrst example was provided by a superconducting source in¬ 
jecting currents into two normal leads OE]. Surprisingly, the effect persists 
even for many-channel conductors [91 HQ]. Furthermore, strong Coulomb in¬ 
teractions can also impose positive cross-correlations in multi-terminal quan¬ 
tum dot systems HBHa. An analysis of the full counting statistics reveals 
a dynamical bunching effect as the origin of the positive cross-correlations 

mm- 

Because of their fascinating properties and potential applications in topo¬ 
logical quantum information [TTl [151 HOl HB HH], Majorana fermions in con¬ 
densed-matter systems have attracted a great deal of interest. However, their 
unambiguous detection in experiments has remained a difficult task: they are 
chargeless and - like the Laughlin quasiparticles in the fractional quantum 
Hall effect - cannot be extracted from their many-body environment. Elab¬ 
orate schemes leading to indirect but conclusive signatures of their presence 
have been proposed and partially realized. 

Recently, several groups uniEQi have reported the identihcation of Ma¬ 
jorana bound states in nanowires by observing a zero-bias peak in tunneling 
spectroscopy experiments. There is, however, no consensus regarding the at¬ 
tribution of this result to the presence of Majorana fermions. The situation 
is similar with the experimental report [21] of a dvr-periodic Josephson effect, 
which cannot yet be unambiguously attributed to the presence of Majorana 
bound states in the Josephson junction. 

There are also many proposals to detect Majorana fermions based on 
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interferometric structures. They can be divided into two classes. The first 
class intends to probe the non-Abelian statistics of Majorana bound states 
trapped in vortices of topological superconductors. In Ref. [221123] the au¬ 
thors study conductance signatures of vortex tunneling in a Fabry-Perot 
interferometer. Another proposal [21| is based on a Mach-Zehnder inter¬ 
ferometer constructed from a topological Josephson junction. In that case, 
Josephson vortices trapping a Majorana bound state propagate along the 
two arms of the interferometer and give rise to a Josephson-vortex current 
A. The presence of absence of another MBS at the center of the interferom¬ 
eter, which can be tuned by a magnetic flux, leads to a striking switching 
between a vanishing A and a nonzero A. The roots of this effect are traced 
back to the non-Abelian exchange statistics of Majorana bound states. 

The goal of the second class of interferometer-based proposals is to hnd 
a signature of Majorana edge states [22l |25l [26l [271 |2H| |29]. They propose 
to use 3D topological insulator heterostructures to build a Mach-Zehnder 
interferometer for Majorana edge states contacted by electronic leads and 
hnd conductance signatures. Reference [26] hnds a signature of chiral Majo¬ 
rana modes in the conductance of a Mach-Zehnder interferometer built in a 
superconductor-quantum spin Hall-superconductor sandwich. 

Biittiker et al. [29] introduced a powerful scattering matrix approach in 
terms of Majorana modes that highlights their special properties and is read¬ 
ily applied to these interferometric structures. 

Motivated by these works, in Ref. [27] this setup was extended to in¬ 
terferometers with two incoming and two outgoing Dirac fermion channels. 
Furthermore it was proposed to study noise correlations in a Hanbury Brown- 
Twiss (HBT) type interferometer, and three signatures of the Majorana na¬ 
ture of the channels were predicted. First, the average charge current in the 
outgoing leads vanishes. Furthermore, an anomalously large shot noise in the 
output ports for a vanishing average current signal is expected. Adding a 
quantum point contact (QPC) to the setup, a surprising absence of partition 
noise was found which can be traced back to the Majorana nature of the 
carriers. 

In view of previous successes of studies of higher-order correlators m 
EHEaE] it is natural to ask [33l |3l] if other Majorana signatures could be 
hidden in higher-order correlations. Therefore we will investigate the full 
counting statistics in multi-terminal structures containing Majorana modes 
to access correlations beyond noise. 

In Ref. [35], the full counting statistics was calculated for a network of 
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Figure 1: Modified Hanbury Brown-Twiss interferometer, see m- Majorana excitations 
will propagate along the boundaries of the two triangular superconducting structures with 
phases cpi, ip 2 - An additional short gapped channel appears at the domain wall between the 
two superconducting regions, forming a quantum point contact characterized by reflection 
and transmission amplitudes r, t. The setup is similar to the one proposed in [25) . 


localized Majorana bound states coupled by tunneling. In this work, in 
contrast, we will use full counting statistics to focus on the transport in 
systems containing one-dimensional, propagating Majorana modes. 

2. Scattering matrix and fnll counting statistics (FCS) 

A very general description of quantum transport in a multi-terminal 
device is to calculate the full counting statistics (FCS) of the transfered 
charges, which contains the full information about the current and higher- 
order current-current correlation functions at zero frequency. Let us denote 
the probability that Ni,N 2 ,... charges are transported into terminal 1, 2,... 
during a fixed measurement time by P{Ni, N 2 ,...). A quantity which is 
equivalent but easier to calculate is the cumulant generating function (CGF), 

dehned by lnx(Ai, A 2 , • • •) = ln(exp(iAiiVi-(-zA 2 iV 2 -I -)), where the average 

denotes the statistical average with the probability P. 

FCS helps to understand transport because it allows an identihcation of 
elementary transport events by decomposing the CGF into a sum of multino¬ 
mial distributions, where the sum indicates the independence of the various 
transport processes in the long-time limit. 

Our goal is to study Majorana interferometers like, e.g., the Hanbury 
Brown-Twiss interferometer with the additional QPC, see Fig. [^and [27] . 

In the absence of interactions, the FCS can be obtained from the 4x4 
scattering matrix S that relates the outgoing Dirac fermions to the incoming 
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Dirac fermions in the following way 


/ '02,e \ / 01,e \ 

02,/i _ ^ 01,/i 

04,e 03,e 

\ 04,/x / V 03,fe / 


( 1 ) 


where 1, 3 label the incoming channels, 2,4 the outgoing channels, e denotes 
an electronic mode, and h denotes a hole mode. 

The matrix elements of the scattering matrix represent the probability 
amplitude of different “processes”. For instance, the matrix element sn is 
the probability amplitude that an incoming electron in lead 1 goes out as 
an electron in lead 2 ; and S 34 is the probability amplitude that an incoming 
hole in lead 3 goes out as an electron in lead 4. The matrix elements already 
contain all the interference effects of single-particle states. 

Two-particle processes: To obtain the probability amplitudes for two- 
particle states, we must take anti-symmetric combinations of products of 
matrix elements. For instance, the probability amplitude that two incoming 
electrons in leads 1 and 3 go out as electrons in leads 2 and 4 is given by 
S 11 S 33 — S 13 S 31 , i.e., the probability amplitude that the incoming electron in 
lead 1 goes out as an electron in lead 2 AND the incoming electron in lead 
3 goes out as an electron in lead 4, MINUS the probability amplitude that 
the incoming electron in lead 3 goes out as an electron in lead 2 AND the 
incoming electron in lead 1 goes out as an electron in lead 4. These two 
processes interfere and lead to effects such as the two-particle Aharonov- 
Bohm effect [36] . 

Three-particle processes: It is also possible to have three-particle pro¬ 
cesses: an incoming lead populated by both an electron and a hole, the other 
incoming lead populated by either an electron or a hole. 

Four-particle process: There is only one such process: all the incoming 
states and all the outgoing states are totally hlled. This (trivial) process 
occurs with probability amplitude det S' = 1 (or probability | det Sp = 1). 

To obtain the full counting statistics we sum up the contributions of all 
these coherent processes weighted by the probability of the input states (oc¬ 
cupation of the incoming leads). For instance, the probability Pij of obtaining 
two particles in the outgoing modes i,j G {1,... 4} is 



k,l=l 
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where k, I run over the incoming modes, and is the occupation of the 
incoming mode k. The indices a and b are the two remaining incoming 
modes (such that the set {a, b, k, 1} = {1, 2, 3,4}). The factor 1/2 is inserted 
to avoid double counting. 

A measurement of the current and its correlation function will not make 
is possible to determine all entries of the matrix because, e.g., processes 
with no charge going out at all are not distinguishable from processes where 
both an electron and a hole go out into the same lead. Therefore, the FCS 
is characterized by only nine independent probabilities: 

• Poo, no net charge in the outgoing leads, 

• P±o, an electron (hole) in the outgoing lead 2, 

• -Poit, an electron (hole) in the outgoing lead 4, 

• P±±, an electron (hole) in both outgoing leads, 

• -P±=F) an electron (hole) in outgoing lead 2 and a hole (electron) in 
outgoing lead 4. 

All of these probabilities can be conveniently retrieved from the scattering 
matrix and the occupation of the leads. The cumulant generating function 
is lnx(A 2 , A 4 ), where A^ is the counting held of outgoing charge in lead i, and 

x(A2,A4)= . (3) 

S2,S4 = -,0,+ 

Since we treat the leads as free fermion reservoirs the knowledge of the prob¬ 
abilities allows to directly access the FCS. The cumulant generating function 
can be alternatively obtained from the Levitov-Lesovik determinant formula 
pUl 137] . which can easily be shown to lead to the same result. 

3. FCS of the Hanbury Brown-Twiss interferometer 

To analyze the full counting statistics of the Hanbury Brown-Twiss inter¬ 
ferometer of Majorana fermions shown in Fig. we insert the appropriate 
expression for the scattering matrix IZ71 


/ ^>2,6 ^ 


( rji — t rji + t r r \ 



V’2,h 

1 

Vi +1 Vi ~ t 



^>4,6 

“ 2 

r -r -r ]2 + t r ]2 + t 


^3,e 

V / 


\ r -r ri 2 + t -r ]2 + t J 


\ ^3,h / 


6 








into Eqs. ([^ and (§. The probabilities are related to charge transfer pro¬ 
cesses of the structure. The probability Po,o does not contribute to the sum 
in Eq. (|^ because it corresponds to processes where no net charge is trans¬ 
mitted to either lead 2 or 4. The remaining probabilities are given by 

F±,o - 20 / 37 ^ ± Re( 7 i)(a - / 3)(475 -M) ( 7 ^ - a/3) , 

Fo,± - 2 a /375 ± Re{r] 2 ){j - 5)(4a/3 + 1) (a/3 - 7 ^) , 

P±,± + a/37<3 - ^(a + /3)(7 + 5) + ^Re( 7 i 72 )(a - /3)(7 - 5) 

± (Re(? 7 i)(a - /3)(1 - 475 ) Re(? 72)(7 - <^)(1 - 4a/3)) 

O 

+ (Re( 7 i 7 ;)(a - /3 )(7 - h) - (a - 7)(/3 - h) - (a - 5)(/3 - 7 )) , 

P±,T + «/^7<3 - ^(a + /3)(7 -F (5) - ^Re( 7 i 72 )(a - /3)(7 - 5) 

± h (Re( 7 i)(a - /3)(1 - 47 ( 5 ) - Re( 72)(7 - <3)(1 - 4a/3)) 

- (Re( 7 i 7 ;)(a - / 3)(7 - 5) + (a - 7)(/3 - (5) + (a - 5)(/3 - 7 )) , 

(5) 

where we introduced the symmetrized occupation numbers a = 4 — me, 
/3 = 4 — nih, 7 = 1“ and <5 = | — ush to simplify the expressions. Since 
the sum of all probabilities has to be equal to one, one hnds 

Fo,o =1 + ^(a -f- /3)(7 -F h) - 1 - AaP'yS 

+ h‘^[{a-'y){/3-6) + {a-6){/3-'y)]. (6) 

The probabilities Ps 2 ,si take into account all the physical processes that 
may occur in the interferometer. Processes with Sj = ±1 correspond to an 
outgoing electron (hole) in lead i, while processes with Sj = 0 correspond to 
processes with no outgoing particles in lead i or to processes with an outgoing 
electron and hole in lead i. 

To better understand the probabilities it is useful to distinguish Eve 
classes of processes depending on the number of incoming particles they in¬ 
volve: the trivial 0 -particle process that contributes solely to Fq^o; 1 -particle 
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processes that contribute to Po,± and P±,o; 2-particle processes that con¬ 
tribute to P±,±, P±,=F, and Po,o; 3-particle processes that contribute to Po,± 
and P±,o; and the trivial 4-particle process that also only contributes to Po,o- 


3.1. Results for the cumulants 

To get all cumulants of the outgoing currents, we have to consider the 
function 


1 

h 


dE Iny 


1 /■“ 

- / dP In 

h Jo 


1 + 


S2,«4 


S2,S4 


ie(s2A2-|-S4A4) 



(7) 


The cumulants C{m, n) are given by the derivatives of this function at A 2,4 = 

0 , 

C{m,n)= {-idx 2 )"'{-idxXT [ dP Iny . (8) 

^ do J A2=A4=0 

As a check we can rederive the results of m For example, the average 
current in lead 2 is given by 


ih) =C(1,0) = - / 


e 

h 

e 

h 

e 

h 


S 2 S 4 

dP (P++ + P+0 + P+- - P-+ - P-0 - P-] 
dE tRe{rii)vi 
dE tRe{T]i){nie - nih). 


(9) 


making contact to the previously obtained result. Similarly the current- 
current cross-correlations read 


^24 


■2 poo 

P(l,l) = ^ / dE ^S2S4Ps2, 

^ S2,S4 

^ ^ 'S2-Ps 2,S4 j ( ^ ^ SaPs2,S4 j 

S2,S4. / \S2,S4 / 


e^R 


dE (ui,e - ni,/j)(u3,e - U3,/i)Re(pi?72) 


( 10 ) 






Figure 2: Cumulants C{m,n) as a function of t (transmission amplitude of the QPC) at 
zero temperature and for 771 = 1, 772 = —1. The vertical axis is in arbitrary units and we 
set Vi = V 3 = e = h = 1. 

in agreement with the results of |27j . 

To go further, let us now specify a set of parameters which capture the 
most interesting physics. Temperature and interferometer asymmetry do not 
bring any new effects, but only smear out certain quantities. We will thus 
set rii ^2 = ±1, and T = 0. Moreover we will adopt a symmetric voltage 
conhguration where Vi = V 3 = V > 0 to access the nontrivial multiparticle 
processes. A plot of the resulting cumulants is shown in Fig. Two obser¬ 
vations can be drawn from this plot. First, it appears that there is a relation 
of the form |C(m,n)| = |C(m-f-n,0)|. Second, the cumulants C(n,0) are 
connected to the cumulants of the binomial distribution as discussed in the 
next section. 

The relation |C'(m,n)| = \C{m + n, 0)| holds because we treat a special 
case : hrst Pq± = P±q = 0 and second, for 7712 = ± 1 , it turns out that either 

P++ = P = 0 or P+_ = P_+ = 0. It then follows that S'(A 2 , A 4 ) is either a 

function of A 2 — A 4 or a function of A 2 + A 4 . The cumulants of the k-th order 
C{m,n), m + n = k, obtained by taking m derivatives of S with respect 
to A 2 and n derivatives with respect to A 4 , are thus exactly the same in the 
latter case (function of A 2 + A 4 ) and only change sign in the former (function 
of A 2 — A 4 ). 

3.2. Half-charge transfers 

The cumulants C(? 7 ,, 0 ) are obtained from the cumulant generating func¬ 
tion 

Inx(A) = In (1 + P+ (P" - 1) + P_ (e-*" - l)) , (11) 
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which follows from the full cumulant generating function (|^ by summing over 
the outcomes in lead 4, and where P± = (1 ± rjitY/A are the probabilities 
of transfering an electron or a hole into lead 2. Note that we still assume 
T = 0, rji^2 = if, Vi = = V >0. 

We may add the side remark that the CGF in the case of f = 0, viz. 
a simple Mach-Zehnder interferometer, corresponds to a trinomial process 
with equal probabilities 1/4 for an electron to electron or to hole transfer 
process. This has been noted in 123. where the corresponding predictions 
for the conductance (it vanishes) and the noise (is quantized with a Fano 
factor 1/4) have been obtained. Recently a detailed analysis of the FCS has 
confirmed this prediction [M] and related it to the topological nature of the 
Majorana mode. 

We are now ready to make the connection with the binomial process. 
Basically, it turns out that 

Inx(A) =21ti (^1 + _ i) + L^(,-iA/2 _ (12) 


This means that we can interpret unit-charge transfer processes as two in¬ 
dependent half-charge transfer processes (notice the 2 multiplying the loga¬ 
rithm, and the factor A/2 in the exponents). In other words, two independent 
binomial processes occuring with probability p = {1 + Vi)/‘2- 

However, we do not expect physical half-charge transport processes be¬ 
cause the existence of charge e/2 quasiparticles in this system is unlikely. We 


can also write Eq. (12) in the equivalent form 


In x(A) =iX -I- 2 In ( 1 -t- 


1 


(e-*^ - 1) 


(13) 


which has a more mundane interpretation. Equation (13) represents two in¬ 


dependent conversion processes: two incoming electrons can either go out as 
electrons with probability y/P+ each, or be converted to a hole with proba¬ 
bility ^/PI. 

Both interpretations, with half-charges or with conversion processes, are 
deceptive. In reality, the two incoming electrons are not independent. How¬ 
ever, the factorization of the process itself, regardless of its interpretation, is 
of special interest. This is a potential signature of Majorana fermions. We 
would like to point out that this factorization is only valid for a symmetric 
interferometer, or at zero bias. This could be potentially attributed to a 
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special property of the fc = 0 Majorana mode (which is self-adjoint) that the 
fc 7 ^ 0 modes do not possess. 

Similar half-charge full counting statistics have been found in other meso¬ 
scopic transport situations, e.g. a double-barrier structure (see Eq. (38) of 
|3H| or Eq. (3.51) of [M]), or a voltage-driven quantum point contact strongly 
coupled to a charge qubit [ 10 ]. 


4. Factorization of transfer processes 

To shed more light on this question, we will now study this factorization 
property for a general scattering matrix. As discussed before, in some cir¬ 
cumstances, it is possible to factorize the charge transfer processes, such that 
the cumulant generating function takes the form 


lnx(A 2 ,A 4 ) = 21nx'(A2, A 4 ) , (14) 

where the “half-CGF” In y' is expressed in terms of half-charge transfers 

^ie{s2X2+Si\i)/2 (15) 


lny'(A2, A4) = In ^ 

S2,S4 = — iOjH" 


P' e" 

S2S4^ 


Note the factor 1/2 multiplying the counting helds in the exponent. 

To see when this factorization is possible, let us compare {x'Y with y: 


(x'Y = 


E 

S2,S4,t2,t4 


P' P' e 


2e(A2 (-52H“i2)/2H-A4(s4H-t4)/2) 


(16) 


The hrst conclusion we draw is that we cannot mix 0-charge transfers with 
1/2-charge transfers in y' in the same outgoing lead. Otherwise, we would 
still have 1/2-charge transfers in (y')^. This means in Eq. (15) we must 
have either 52,4 = ±1 or § 2,4 = 0). The case of 0-charge transfer only is 
not interesting, so we focus on purely half-charge transfers. We obtain the 
following equations for the probabilities 

{P'^^f = , 

= P±T , 

‘^P±+P±- = P±o , 

‘^P+±P-± = Po± , 




+ 2P'P' = P 


00 • 


(17) 
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The first two equations fully fix the four probabilities P', the three remaining 
equations are thus nontrivial conditions the probabilities P must satisfy to 
make the factorization possible. 

If we specialize to two-incoming-particle processes, we can rehne the con¬ 
ditions on P. First, if we only have two-incoming-particle processes (this can 
be tuned by choosing proper occupations of the incoming leads, such as a 
symmetric bias and zero temperature), then P±q = Po,± = 0. The half-charge 
transfer process (i.e., the factorization) is then possible if either 

y/P++ + \/ P — = 1 or y/p)!r -l- \/ P-+ = 1 . (18) 


Thus, we are really looking for a property of the scattering matrix (that is 
independent of the incoming leads’ occupation) rather than a property of the 
full counting statistics themselves. Typically, if two-particle processes and 
one-particle processes mix, we have no chance of hnding such a factorization 
of the FCS. 

We now specialize to the case where both electron incoming modes are 
fully occupied, and the hole modes fully unoccupied. This leads to the prob¬ 
abilities 


V P++ = l-SllSsS — SisSsil , 

\/ P — = |s2lS43 ~ ■S23S41I 5 

\/ P+- = |■SllS43 ~ 5135411 , 

V P-+ = |52i533 — S 23 S 31 I . (19) 


Thus, checking the factorization property given in Eq. (18) has been re¬ 
duced to a condition on the scattering matrix. The strategy will now be to 
investigate interferometers of increasing complexity by writing down their 
scattering matrices and checking whether Eq. (18) is fulfilled, i.e., whether 
the factorization property holds. 


4 . 1 . Double point-contact geometry 

We hrst consider the double point-contact interferometer shown in Fig.[^ 
where the labels for the incoming and outgoing Majorana modes are dehned. 
The conversion of Dirac to Majorana modes is achieved by 


/ riA riA 0 0 \ 

1 i —i 0 0 

^2 0 0 -r]c -Vc 

y 0 0 i / 


( 20 ) 
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Figure 3: (a) Double point-contact geometry, (b) Central island geometry, (c) General 

scattering matrix. The numbers 1, 2, 3,4 in the top-left and botton-right corners label the 
four incoming chiral Majorana modes, whereas the numbers in the top-right and bottom- 
left corners label the four outgoing Majorana modes. Incoming and outgoing Majorana 
modes are related by the scattering matrix Sm- 


see 


Appendix A Here y]A-,^c = =tl depend on the parity of the number 


of vortices. The minus sign in front of 7]c takes care of the vr Berry phase 
obtained when a Majorana fermion makes a 2tt rotation. The back-conversion 
is similar, 


’6c 



( -i 

Vb 

0 

0 \ 

1 

i 

Vb 

0 

0 

71 

0 

0 

—i 

Vd 


V 0 

0 

i 

Vd J 


( 21 ) 


The tunneling of Majorana fermions is described by the S'0(4) matrix 


/ 0 


Sm = 


r 

0 


—t 0 r \ 
0 t' 0 
r 0 t 


V 


-t' 0 


( 22 ) 


0 / 


The total scattering matrix is given by the product of these three scat¬ 
tering matrices 


S = S,,SmSc 


(23) 
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which is real. Using Eq. (19), one obtains 



(24) 


(25) 


i.e., the sums of the probability amplitudes are seen to fulhll Eq. (18). Thus, 
factorization holds, and the two-incoming-electron processes can be written 
as two equal and independent processes. 


Central Majorana island geometry 
We define the matrix 


/ 0 0 0 1 \ 
10 0 0 
0 10 0 
\ 0 0 1 0 / 


(26) 


as well as i? = diag(ri,r 2 , rs, r 4 ), and T = diag(fi, ^ 2 , ^ 3 , ^ 4 )- Using (XR)* = 
rir 2 r 3 r 4 = r is proportional to the identity matrix, we find 


Sm = XR -^XT(1 + XR+ iXRf + (XRf)XT . (27) 

1 — r 
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The relevant matrix elements of the total scattering matrix S = 
read 


-Sll 

-533 

Sl3 

-S 3 I 

S 2 I 

"S 43 

"S 23 

S41 


n - 

rs - 

Ti - 

r 2 - 

n - 

r3 - 


-Ta 


-r 2 


tir2r^rAti + t4’^i^2 - - ^ 1 ^ 2 ) 

1 — r 

hrArir2t3 + t2rztA - - htj) 

1 — r 

t4^i^2^3^4 + tir2h - i{t4rir2t3 - tir2r3t4) \ 

1 - r / ’ 

t2r3rArit2 + hvAti - f(t2r3^4ti - t3rArit2) \ 
1-r ) ’ 

tir2r3r4ti - tArit2 + i{tAti + tit2) 

1 — r 

t3r4rir2h - t2r3U + ^(^ 2^3 + ^ 3 ^ 4 ) 

1 — r 

-Urir2r3tA + tir2t3 + i{tArir2t3 + tir2r3t4)\ 


1 — r 


) ' 


-t2r3r4rit2 + hrjti + i(t2^3^4ti + t3r4rit2) \ 
1 — r y 


SbcSMSc 


(28) 


If only one QPC has a non-vanishing transmission, the scattering matrix 
reduces to the HBT case. For t 2 = t, = 1, and ti = = 0, we get back 

the setup HBT with one QPC. For the case where all QPC are equal, t* = f, 
the factorization also holds. The same is true for ti = = t and t 2 = = 0 

although this situation is less symmetric. 

We were unable to prove the factorization property for general values 
ti analytically. We therefore used a numerical approach: in the spirit of a 
Monte-Carlo calculation, we generated many scattering matrices by using 
random values of the transmission amplitudes f*. The factorization was then 

conhrmed by checking that \^P++ + VP _= |'Sii<S 33 — S 13 S 31 I -|- |s 2 i<S 43 — 

-S 23 -S 41 I =0 or 1. This is indeed the case, i.e., the factorization property 
holds for the island geometry as well. 


4-3. Most general scattering matrix 

We now consider the most general case: a general S'0(4) matrix accom¬ 
panied by the appropriate Dirac to Majorana conversion / back-conversion. 
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see Fig. The conversion from 



V 

the back-conversion is similar 

c _ 1 

The total scattering matrix is 


Dirac to Majorana channels reads 


1 

1 

0 

0 

\ 

i 

—i 

0 

0 


0 

0 

-1 

-1 


0 

0 

i 

—i 

/ 


( -i 

1 

0 

0 \ 

i 

1 

0 

0 

0 

0 

—i 

1 

1 0 

0 

i 

1 / 




SbcSnSc , 


(29) 


(30) 


(31) 


where Sm £ 50(4). Again, we demonstrated the factorization property 
numerically by generating a set of random 50(4) matrices and constructing 
the corresponding scattering matrix 5 dehned by Eq. (31). The factorization 
property (18) was found to hold even in this most general case. 


4.4- More than four terminals 

Finally, let us discuss how the peculiar properties of the CGF can be 
generalized to setups with more than four terminals. In the case of six 
terminals, for instance, a general scattering matrix between incoming and 
outgoing Dirac fermions can be expressed as 


5 = SbcSnS ,, (32) 

where Sm is a general 50(6) matrix, and Sc and She denote the 6 x 6 con¬ 
version and back-conversion matrices. If particles enter via the leads 1,3,5 
and exit via the leads 2,4,6, then the most general CGF is a generalization 
ofEq. 

x(A 2 , A4, Ae) = _ ( 33 ) 

■ 525 - 54,56 = —, 0 ,+ 

If we focus again on the limit of equal positive bias voltages lA = V 3 = 
V 5 > 0 , then all incoming electron modes are occupied and all incoming 
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hole modes are empty. In that case, in analogy to Eq. (19), one hnds, e.g., 
for the probability P+++ of Ending particles in all ontgoing modes and the 
probability P++- to hnd particles exiting in leads 2 and 4 and a hole exiting 
via lead 6 , 


■\/ P+++ — ^ |( — 

r&Si35 

\/ P++- = ^ |( —l)*®^^^^Sl,-PlS3,P3S5,P6|- (34) 

'P&Sise 

where the snmmation is over the six permutations of the respective set. With 
our choice of phases in the conversion matrices, one hnds that for an arbitrary 
SO (6) matrix Sm only processes with outgoing charge 3 and —1 have a 
nonzero probability. 


P+++ + P +— + P-+- + P —+ — 1, 

F++_ = F+_+ = F_++ = F___ = 0. (35) 


Using the hrst line to eliminate F+++, one hnds that the CGF is given by a 


generalization of Eq. (13), 


lnx(A 2 , A 4 , Ag) — i{X 2 + A 4 + Ag) + In 


1 + P+—{e 


— 22(A4 + A6) _ 


1 ) 


^p_^_(g-2i(A2+A6) _ 1) + p__^(g-2i(A2+A4) _ 


(36) 


Evidently, for a general scattering matrix the six-terminal result contains 
more free parameters than the four-terminal case. This means that it is no 
longer possible to factorize the ECS in general. 

However, one still hnds that only a very limited number of scattering 
processes are possible: either two of the incoming eletrons are transmitted as 
holes (with probabilities F^_, F__|__ and F_^) or all electrons are trans¬ 

mitted as electrons. 


5. Discussion/Conclusion 

To summarize, we have calculated and analyzed the full counting statis¬ 
tics of a modihed Hanbury Brown-Twiss interferometer for chiral Majorana 
fermions which contain information about higher-order current correlation 
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and generalize the resnlts presented in [27j. Most of the calcnlations in this 
paper are valid only at zero energy, at which particle-hole symmetry enforces 
important constraints on the scattering matrices. Some of the obtained re¬ 
snlts are expected to remain valid at hnite, bnt low, energies. However, we 
also do not keep track of dynamical phases which are important even at low 
energies. 

The fnll connting statistics calcnlated in this paper exhibit an interest¬ 
ing factorization property that points towards the interpretation of a nnit- 
charge transfer process as two independent half-charge transfer processes. 
We checked the factorization property for increasingly more general scatter¬ 
ing matrices and conhrmed nnmerically that it holds for the most general 
S'0(4) matrix at specihc conhgnrations of voltage bias and at zero tempera- 
tnre. It is tantalizing to interpret this property as a signatnre of Majorana 
fermions, however, since we do not expect physical 1/2 charges to be trans¬ 
ferred, its interpretation remains open at present. 
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Appendix A. Majorana interferometers basics 

In this appendix, we list some of the basic bnilding blocks needed to 
stndy concrete interferometric strnctnres, compnte their scattering matrix, 
and check whether the FGS factorization holds. 

Appendix A.l. Dirac to Majorana converter 

This is the bnilding block to connect Majorana fermions to external Dirac 
fermions [221E5]. The scattering matrix that relates chiral Dirac fermions to 
chiral Majorana fermions at T-jnnctions is fnlly determined by particle-hole 
symmetry and nnitarity. Its form is given by 
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such that 


i j - I 

The only remaining freedom is a physically unimportant phase factor that 
corresponds to gauge transformations 

S-conv ^ S-conv f q j • (A.3) 

As long as the Dirac fermions in the external channels do not interfere with 
one another, these factors are irrelevant. 

The back conversion of Majorana to Dirac fermions is described by the 
time-reversed scattering matrix 


—ronv /— 


1 i 
1 -i 


(A.4) 


Appendix A.2. Majorana fermion point contact 

Because of particle-hole symmetry, the scattering matrix that relates chi¬ 
ral Majorana fermions before and after a quantum point contact (QPC), or 
any form of tunneling, must be real. The scattering matrix is therefore a ro¬ 
tation matrix that can be parametrized by a single angle. We use r, t as the 
(real) reflection and transmission amplitudes of the point contacts to write 


r t 
—t r 


(A.5) 


Note that tunneling between a pair of Majorana channels directly after 
the conversion from Dirac fermions can be described by the product of two 
scattering matrices 


*S*tot Atunnel *S*conv 


r + it r — it 
—t + ir —t — ir 


y ie*" —ie 

where we dehned r -(- ft = e*". We can write this as 


A+ot — Snj 


(A.6) 


(A.7) 


The effect of tunneling is therefore equivalent to a redehnition of the phases 
of the incoming Dirac modes and can be disregarded. 
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